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SUMMARY 


In  this  report,  the  problem  of  reconstructing  the  power  spectral  density  of  a 
Gaussian  signal  process  from  hardlimited  observations  taken  at  Poisson  sampling  instants 
was  considered.  Hardlimiting  the  signal  represents  the  most  drastic  form  of  amplitude 
quantization  that  is  possible  - but  yet  allows  meaningful  analysis  to  be  performed.  The 
sampled  hardlimited  data  is  represented  with  only  one  bit  of  information  which  gives  a 
significant  reduction  in  the  amount  of  storage  area  needed  for  retaining  the  data.  Also,  the 
processing  of  the  data  becomes  simplified  because  of  the  simple  representation.  By 
Theorem  1 , the  power  spectral  density  estimate,  is  shown  to  be  asymptotically  un- 

biased as  the  number  of  observation  points  approaches  infinity.  The  asymptotic  rate  of 
convergence  for  the  bias  (of  the  estimate)  is  identical  to  the  rates  of  convergence  for  the 
cases  of  Poisson  spectral  estimation  (without  hardlimiting)  and  of  periodic  spectral  estima- 
tion (to  the  aliased  spectral  density).  Furthermore,  because  the  estimate  is  asymptotically 
unbiased,  no  aliasing  of  the  spectral  density  occurs  in  spite  of  the  hardlimiting. 

The  covariance  of  the  estimate  was  shown,  however,  only  to  be  bounded  as  the 
number  of  observations  tends  to  infinity.  Although  not  indicated  by  the  analysis  up  to  this 
point  (primarily  Theorem  1 ),  the  number  of  hardlimited  samples  needed  for  the  same 
quality  of  estimate  of  the  spectral  density  may  increase  significantly  over  the  number  of 
samples  needed  in  the  case  when  no  hardlimiting  is  done.  The  focus  of  future  work  will  be 
to  obtain  tighter  bounds  on  the  covariance  of  the  estimate  and  to  show  that  these  bounds 
tend  to  zero  as  the  number  of  observations  tends  to  infinity. 


Q O 


INTRODUCTION 

This  paper  is  concerned  with  the  estimation  of  the  spectral  density  function  of  a 
continuous-time  stochastic  process  from  observations  taken  at  discrete  instants  of  time  after 
the  process  has  been  hardlimited.  In  [ 1 ] , a continuous-tinje  estimate  using  the  hardlimited 
realization  of  a stationary,  Gaussian  stochastic  process  was  shown  to  be  a consistent  estimate 
of  the  spectral  density  function  of  the  original  process.  Hinich  [2]  deals  with  the  consistent 
estimation  of  aliased  versions  of  nonbandlimited  spectra  from  hardlimited  Gaussian  processes. 
Brillinger  [3]  provides  consistent  estimates  of  the  cross-spectrum  and  spectra  of  stationary 
bivariate  Gaussian  processes  from  the  zeros  of  the  processes.  Finally,  Masry  and  Lui  [4]  and 
Masry  [5]  provide  consistent  estimates  of  spectra  of  continuous-time  processes  from  obser- 
vations taken  at  Poisson  sampling  instants. 

In  this  paper  we  consider  the  problem  of  estimating  the  spectral  density  function  of 
a Gaussian  process  - the  estimate  is  based  upon  hardlimited  (+1,  -1  or  0)  samples  of  the 
process  taken  at  Poisson  sampling  instants.  The  proposed  estimate  is  shown  to  be  asymp- 
totically unbiased.  For  the  covariance  of  the  estimate,  all  but  two  terms  are  shown  to  go  to 
zero  and  these  terms  are  easily  shown  to  be  bounded.  Future  work  includes  finding  bounds 
for  the  rate  at  which  the  bias  asymptotically  goes  to  zero,  along  with  showing  that  the 
covariance  goes  to  zero,  and  the  rate  at  which  the  covariance  tends  to  zero. 

Let  X =|X(t),  - °°<t<<»}bea  real  fourth-order,  stationary  Gaussian  process  with 
zero  mean,  correlation  function  R(t)  and  spectral  density  function  $(X).  Define  the  hard- 
limited  version  of  the  process  X(t)  as  Y(t)  = sgn  X(t)  where  sgn  X is  equal  to  1,-1  or  0 
according  to  X being  positive,  negative  or  zero,  respectively.  Then  Y(t)  is  a fourth-order,  zero 
mean,  wide-sense  stationary  process  with  correlation  function  [3,6 J 

RY(T)  = E[Y(t  + T)  Y(t)l 

= (2/jr)  Arcsin  (R(r)/R(0)).  ( 1 ) 

Since  all  amplitude  information  is  lost  when  the  X process  is  hardlimited  we  assume  that 
R(r)  is  normalized  such  that  R(0)  = 1 . 

The  following  relationships  hold: 


R(t)  = 


/ 


Ry(t)=  j <P y(X)eitXdX 


and  if  R(t)  eLj  and  Ry(t)  eLj  then  the  inverse  relationships  of  (2)  hold 
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T^" 


oo 

0(X)  = / R(t)  e”itX  dt/2w 

— OO 

OO 

0Y(X)=  / Ry(t)e-itX  dt/2?r . (3) 

— OO 


Using  (1)  and  (3)  and  noting  that  Ry(t)  is  an  even  function,  the  spectral  density  function  can 
be  written  as 

OO 

0(X)=  J sin  [it  Ry(t)/2]  cos  tX  dt/w.  (4) 

o 

The  problem  that  is  considered  here  is  to  estimate  0(X)  by  sampling  the  hardlimited  process 
Y(t)  at  discrete  instants  of  time. 


which  are  determined  by  a Poisson  point  process. 

The  proposed  Poisson  finite  sample  estimate  of  0(X),  using  K'  discrete  observations  of 
Y(t)  which  are  taken  from  a Poisson  point  process,  is 


0hj(X)-ff«  ^ wN(tn)  cos  Xtn  sin  2n  ^ ^lk+n^^k^ 
n=l  k=l 

The  |tn|  are  generated  by  a Poisson  point  process  on  (0,<*»),  i.e., 


(5) 


tD  = o 


tn  = tn-l+0;n  n=  1,2,3,...  (6a) 

where  the  |an|  are  independent,  identically-distributed  (positive)  random  variables  with  a 
common  exponential  distribution  F(x)  = 1 - e~0x  and  P is  the  average  sampling  rate.  The 
process  T = |tn|  is  assumed  to  be  independent  of  the  process  X(t).  The  probability  density 
function  of  tk+n  -t^  is  independent  of  k and  is  given  by 


fn(t)  = P 


(n  -1 )! 


e-^,  t>0. 


(6b) 
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The  covariance  averaging  kernel,  wj^t),  is  given  by 


where 


OO 

h(t)  = f H(X)  eitX  dX 

— OO 


and  H(X)  eLj  is  even  and  integrates  to  1.  The  spectral  window,  Wn(X),  is  defined  by 
Wn(X)  = Mn  H(MnX) 

where  < N is  such  that  Mjq  -*•  °°  and  Mjq/N  -*•  0 as  N -*• «>. 


The  primary  results  of  this  paper  are  concerned  with  the  asymptotic  properties  of 
the  spectral  estimate  (5).  Theorem  1 establishes  that  the  estimate  is  asymptotically  un- 
biased as  the  number  of  observations  N tends  to  infinity.  Furthermore,  under  mild  con- 
ditions of  integrability  on  the  correlation  function  (tR(t)  eLj),  the  asympototic  rate  of 
convergence  of  the  estimate  to  the  true  spectral  density  (of  the  original  nonbandlimited) 
signal  process  is  0(1 /N),  i.e.,  the  bias  error  tends  to  zero  at  the  rate  of  at  least  1/N.  As  a 
consequence  of  Theorem  1,  the  estimate  ^jq(X)  eliminates  “aliasing”  of  the  spectral  density 
(the  folding  of  high  frequency  power  to  lower  frequencies  which  occurs  with  periodic  samp- 
ling) in  defiance  to  the  hardlimiting.  We  also  note  that  the  asymptotic  rate  of  decrease  of 
the  bias  is  identical  to  the  asymptotic  rate  of  decrease  of  the  bias  for  Poisson  spectral  esti- 
mation when  no  hardlimiting  occurs,  which  is,  in  turn,  identical  to  the  rate  of  decrease  of 
the  bias  of  continuous  time  spectral  estimation  [5] . The  covariance  of  the  estimate  is 
shown,  however,  only  to  be  bounded  as  the  number  of  observations  tends  asymptotically  to 
infinity  .The  bounding  of  the  convariance  is  done  in  a series  of  lemmata,  Lemma  2 through 
Lemma  5.  We  show  that  the  covariance  is  (asymptotically)  bounded  by 


Cov 


R2(u)  du 
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+ ( J I Ry(u)  I du  ) | + 0(A)0(oj) 


We  note  that  since  the  covariance  of  the  estimate  - at  least  by  the  analysis  contained  in  this 
report  - does  not  tend  to  zero,  but  is  only  bounded,  the  question  as  to  whether  or  not  the 
estimate  is  consistent  is  still  open. 


THE  JOINT  DENSITY  FUNCTION  FOR  POISSON  SAMPLES 

Before  stating  the  main  results,  we  collect  together  in  this  section  some  preliminary 
results  which  will  be  helpful  in  the  derivation  of  the  bias  and  variance  of  the  estimate  (5). 
We  begin  by  defining  Ry(N)(n)  as 


1 1 — it 

RY(N)(n)=  (1/N)  J Y(tk+n)Y(tk) 


and  the  estimate  (5)  becomes 


?N(X) = ^ X wN(tn>  cos  Xtn  sin  ’rRy(N)(n)/2  • 
n=l 

Note  that  Ry^(n)  is  not  an  estimate  of  Ry(n)  since 

OO 

E Ry^N)(n) = (l-jj)  / Ry(t)fn(t)dt 


and  as  N becomes 


+ 0) 
J2 


max  (1,  b(0)] 


(/ 


b(u)  du 


oo 

/ RY(t)fn(t)dt  , 
o 


which  in  general  is  not  equal  to  Ry(n)  but  can  be  considered  as  an  estimate  of  Ry(tn). 

The  following  identities  are  easily  established  (recalling  that  R(0)  = 1)  from  (1) 

sin  [7r Ry(t)/2]  = R(t) 

o 1/2 

cos[7rRy(t)/2]  =[1  -R2(t)]  . (9) 

These  identities  are  used  later  and  also  here  to  simplify  the  two  term  Taylor  series  expansion 
of  sin  [it  Ry(N)(n)/2]  about  the  point  7rRy(tn)/2 

sin  [7T  Ry^(n)/2] 


= R(tn)  + An  jt  | Ry<N>(n)  - Ry(tn) 


1 - Rz(tn) 


1/2 


12 


(10) 


where  |An|  < 1.  The  Taylor  series  expansion  is  used  to  determine  the  asymptotic  behavior 
of  the  covariance  of  the  spectral  density  estimate. 

Next  we  borrow  a lemma  from  [4] . 


Lemma  1:  Consider  the  integrals 

OO. 

dn  = J g(x)  fn(x)  dx  , n = 1,  2, . . . 

o 

where  g eL  j and  fn  is  given  by  (6b).  Then 

oo  oo 

(a)  dn=0  J g(x)dx<~ 

n=l  o 

and,  if  in  addition  x g(x)  e Lj , then 

OO  OO 

(b)  ^ndn  = 0 j (1  +0x)g(x)dx<«>  . 

n=l  o 

The  proof  is  straightforward  and  given  in  [4) . 
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The  result  of  the  following  proposition,  the  joint  density  function  of  the  sampling 
points,  is  used  in  the  derivation  of  the  covariance  of  the  estimates  (5)  and  is  independently 
of  interest. 


Proposition.  Let  |tn}  be 
fn,n+k(x’y)> is  given  by 

*n,n+k(x,y^  “ 
for  n = 1,  2, . . .,  and  k = 


defined  by  (6a).  Then  the  joint  density  function  of  tn  and  tn+k, 

fk(y  - x)  fn(x)  y > x 

0 y 

1,2,...,  where  fn(x)  is  the  density  function  (6b)  of  tn. 


Proof.  The  following  two  relationships  are  used  in  the  derivation  [6] 

*n,n+k(x,y)=  ^n+k^^n  “ x^n^x^ 

fn+k(y|tn=x*  = ify  Fn+k(y|tn~x) 

where  Fn+k(y|tn  = x)  is  the  conditional  distribution  function  of  tn+k  given  that  tn  = x. 
From  equation  (6a),  tn+k  is  given  by 


*n+k  ln  + ^n,k 


where 


n+k-1 

An,k=  S “i  • 
i=n 

Now,  by  recalling  the  fact  that  {crj } are  independently,  identically  distributed 

Fn+k(y ltn  = x>  = P^n  + \i,k  ^ y'  *n  = x> 


= p(An,k<y-x) 

= Fk(y  - x)  . 

Using  (1  lb),  fn+k(yltn  = x)  = fk(y  - x)  and  substituting  this  into  (11a)  gives  the  desired  result. 
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BIAS  OF  THE  ESTIMATOR 


In  this  section  we  show  that  the  bias  of  the  estimator  goes  to  zero  asymptotically  as 
N -*•  <».  The  bias  of  the  estimator  (5)  is  defined  as 


b[^jsj(A)]  - E [<£n(A.)|  - 0(A)  . 


The  following  theorem  establishes  that  the  estimator  is  asymptotically  unbiased. 


Theorem  1.  Let  R(t)  e Lj.  The  estimate  is  an  asymptotically  unbiased  estimate  of 
0(A),  i.e., 

(a)  b[$N(A)]  =o(l) 

uniformly  in  A as  N -*•  Furthermore,  if  t R(t)  e L j , then 

(b)  b [?n(A)|  = 0(1/N)  , 


uniformly  in  A as  N -*• 


Proof.  By  (10)  the  estimator  becomes  equal  to 


where 


?NW  = Ql(X>N>  + Q2<X’N> 


Qj(A,N)  = ^ ^ wN(tn)cos  Atn  R(tn) 


2„  ,1  l'2 


Q2(A,N)  = -^-  AnwN(tn)  cos  A tn  RY(N)(n)  - Ry(tn)  1 - R2(tn) 


Now  taking  the  expected  value  of  Qj(X,N)  (with  respect  to  T = {tn}  since  the  Y process 
does  not  appear) 

oo  N 


EQj(A,N)  = ^ J wN(t)  cos  At  R(t)  ^ fn(t)  dt 


I 


and 


N 


lim  wN(t)  S fn(t)  = 0 
N^°°  n=l 

since  lim  wj^(t)  = h(o)  = 1 and  by  Lemma  1 lim  Z fn(t)  = 0.  Also  note  that 


N 

|wN(t)cosXt  ^ fn(t)|<0  -Ha 
n=l 

where  Ha  = / 1 H(A)  I dX.  Therefore,  by  the  Lebesgue  dominated  convergence  theorem 

OO  N 

lim  EQj(X,N)  = ^ f lim  wN(t)  V fn(t)  cos  Xt  R(t)  dt 
XI  ^ ^ o N^°°  n = 1 


N 


= — f cos  Xt  R(t)  dt  = 0(X)  . 
7 r j 


(14) 


To  take  the  expected  value  of  Q2(X,N),  the  expectation  is  first  taken  with  respect  to  the 
Y process  and  then  with  respect  to  the  T process.  Taking  the  expected  value  with  respect 
to  Y yields  the  bound 


N N-n 

IEYQ2(X,N)l<Ha  ]T  jj;  £ RyOktn-tk) 
n = l 


k=l 


- MW 


(15) 


since  |wN(t)|  < Ha  and  [ 1 - R^(t)]  ^ < 1.  Now  taking  the  expected  value  of  (15)  with 
respect  to  T gives 


|EQ2(X,N)|< 


*|t 

•n=l 


(l'§)  / Ry(t)  fn<t)  dt  - / RY(t)fn(t)dt 


H 


_a 

N 


N <» 

£ n f Ry(t)fn(t)dt 


(16) 


n=l  o 
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and  by  the  Kronecker  lemma  [6] , since 


N 00 


Urn  ^ f RY(t)fn(t)dt  = 0 / Ry(t)  dt  < 00 , 


the  expected  value  of  the  Q2(X,N)  term  goes  to  zero: 

iEQ2(X,N)|-»OasN^oo  . 

Therefore,  the  estimator  (5)  is  asymptotically  unbiased  and  part  (a)  is  established. 
To  establish  part  (b),  the  bias  can  be  written,  using  (12)  through  (14),  as 


b|0N(X)]  =b,(X,N)  + EQ2(X,N) 


where  bj(X,N)  is  defined  by 


bj(X,N)  = -^  J wN(t)  cos  Xt  R(t)  ^ fn(t)  dt  . 


n=N+l 


en=  f |R(t)|fn(t)dt 


n=  1,2, ... 


(nen)' 


is  summable  by  Lemma  1 (b).  Hence, 


lb,(W = 


and  by  (16) 


|E  Q2(X,N)|  = 0(1 /N)  . 


Therefore, 


b|?N(X)]  =0(1/N) 


and  part  (b)  is  established. 


Q.E.D. 
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COVARIANCE  OF  THE  ESTIMATOR 


The  covariance  of  the  estimator  (5)  is  given  by 

Cov  j ^(X),  0n(w)J  = E | ^f(X)  - E$jyj(X)J  - E^"N(o;)j 

By  replacing  E^n(X)  with  the  bias  expression,  along  with  noting  that  the  bias  is  not  a 
random  variable,  the  covariance  becomes 


CovI^XX^gj)!  ~ E [$n(X) - 0(X)|  |^(w)-0(oj)j-b|?N(X)j  bJ^N(w)j  . 

Since  the  bias  goes  to  zero  asymptotically  as  N -*■  <»  we  need  only  to  investigate  the 
asymptotic  behavior  of  E (?n(X)  - 0(X)1  [0n(o>)  - 0(«)1 . By  using  the  two  term  Taylor 
series  expansion  (10)  of  sin[jrRY^(n)/2]  we  have 


:|^n(X)  - m\  - 0(oj) 


= E 


N 


fa  ^ wN(tn)  cos  X tn  R(tn) 

n = l 


N 


Jp  ]>  AN  cos  Xtn[  Ry^V)  - RY(tn)]  [ 1 - R2(tn)] 1 /2  - 0(X) 


n=l 


N 


fa  ^ wN(tn)coscotnR(tn) 
n = l 


N 


Jp  BN  wN(tn)  cos  cotn  R^OO-RyCt,,)  l-R2^)  1/2-tf>(w) 


where 


TjCN.X.w) 


Ei 


N 

fa  ^ WN(tn>  cos  Xtn  R(‘n)  ~ «X) 
n=l 


N 

fa  2 wN<tn>  cos  wtn  R^n>  ~ 
n=l  j 


T2(N,X,w)  = E 

N 

20  J ANwN<*n)  cos  Xtn 

Ry^(n)  - RY(tn) 

1-R2(tn) 

1/2 

n=l 

■ 

- 

! N 

i v-1 

~fa  2 wN^n)  cos  wtn  Wtn)  “ #(«) 
n = l 


T3(N,X,w) 


N 

1/2 

E 

Jp  ^ BNwN(tn)cos<otn 

R^)(n)-Ry(tn) 

1 — R2(tn) 

n=l 

• 

- 

! 

I 


I 


N 


fa  ^tn)008^^)-^) 
n = l 


T4(N,X,co)  = E 


N 

2p  1 ANwN^n>  cos  Xtn 

R™  (n)-Ry(tn) 

1-R2(tn) 

1/2 

n=l 

. 

■ 

N 


JP  ^ BN  wN^  cos  wtn 


n = l 


R^N)(n)-RY(tn) 


1 -R2(tn) 


1/2 


The  asymptotic  behavior  of  each  Tj(N,X,  w)  will  be  obtained  in  a series  of  lemmata.  We 
begin  with  the  T2(N,X,«)  and  T3(N,X,w)  terms. 


Lemma  2.  Let  R(t)  eLj  and  Ry(t)  eLj.  Then  T2(N,X,to)  goes  to  zero  asymptotically  as 
N goes  to  infinity. 

Proof.  The  second  term  of  the  product  in  T2(N,X,cu)  is  (recall  (12))  equal  to 
««)  Q2(X.N) 
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and  by  ( 17)  goes  to  zero  asymptotically  as  N ->  «>.  The  first  term  of  the  product  is,  after 
taking  the  expectation,  equal  to 


'N  f f 

fa  J J WN<S>  wN(t> cos  cos  Ry*®)  f 1 - r2(*)1  R(t) 


o o 


N N 

’2  I 5 fn,k(s>t)dsdt 
n=l  k=l 

where  fn,k(s>t)  is  8iven  by  the  Proposition.  This  term  is  bounded  by 

^ oo  oo  N N 

f f RY(S)  R(t)  ± V £ n fn>k(s,t)  ds  dt  . 
° ° ‘ n=l  k=l 


Consider  the  sum 


N N 


’S,t^  2 2 fn,k(s’1 


n=l  k=l 


N-l  N-n 


2 2 <f».  n+k(s,t)  + fn  n+k(t,s)J  + V fn(min|s,t|) 

« < « 


n=l  k=l 


for  the  case  s<t.  Then 


N-l  N-n 


S(N,s,t)  = ^ £ fk(t-s)  fn(s)  + 2 fn(s) 


n=l  k=l 


< 2 f°(S)  2 fk(t_S)+  2 fn(S> 


which  converges  to  0(1  +0)  as  N -*• <*>  and  by  the  Kronecker  lemma  [61 


N N 

n=l  k=l 


The  same  argument  holds  for  t < s and,  therefore,  by  the  Lebesgue  dominated  convergence 
theorem 

oo  oo  N N 

/ /RYtoWOjL  2 2 n fn>jc(s,t)  ds  dt  -*•  0 

° ° n=l  k=l 

as  N -*•  oo.  Hence 

T2(N,X,u>)  -*■  0 as  N -*■  °°  . Q.E.D. 

Lemma  3.  Let  R(t)  eLj  and  Ry(t)  eLj.  Then 
T3(N,X,g>)  -*■  0 
as  N -+oO' 


The  proof  of  Lemma  3 is  identical  to  the  proof  of  Lemma  2. 
By  multiplying  the  factors  of  Tj(N,X,oj)  we  obtain 

4 

T j (N,  X,  w)  = ^TU(N,X,«) 
i = l 


where 


Tj  ,(N,X,o>)  = — 1—x  E 

1,1  m )2 


N N 

^ wN(tn> cos  Xtn  R(tn^  * ^ wN(tn)  cos  wtn  R(tn) 
n=l  n=l 


Tl,2(N.X,U)-^lE 


N 


wN(tn)cos  GJtn  R(tn) 

n=l 


T1>3(N,X,<o)  = T1  2(N,w,X) 


Tj  4(N,X,w)  = ^(X)0(o;)  . 
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We  note  that  Tj  ^(N.X.w)  = -0(X)  • Qj(cj,N)  and  by  the  proof  of  Theorem  1 


■ : 


Y 


\L 


lim  T j 2(N,X.co)  = -^(X)  0(w) 
N ->•  °o  ’ 


(18) 


The  term  T j > j (N, X, go)  is  divided  into  three  sums: 


T,  ,(N,X,w)=  ^Tj  j^N.X.oj) 
i=l 


where 


N 


t11,(N,X,w)  = ^2  E £ wj,(tn)  cos  Xtn  cos  wtn  R2^) 


n = l 
N 


Tl,l,2<N’X’w>=-“2  E J £ WNa^WNa^cosXtnCOStotkRfVR^) 


n = l k>n 
N 


t1,1,3(N’X’w)='^2E  J S wN<tn>wN<tk)co,Xtnco,wtkR(tn)R(tk) 


n=l  k<n 


The  Tj  i j term  will  be  disposed  of  first.  By  taking  the  expected  value  (with  respect  only 
to  the  j tk  | ) 


N 00 

Tl,l,l(N»x»w)  - — j ^ J wj,(t)  cos  Xt  cos  cut  R2(t)  fk(t)  dt 


k=l  0 

and  by  the  Lebesgue  dominated  convergence  theorem 


lim  Tl,l,l(N,X,w) J cos  Xt  cos  cut  R2(t)  dt 
N-»°°  0*  o 


(19) 


We  note  that  Tj  j 2 an(l  ^ 1 1 3 can  be  rewritten  as 
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N-l  N-n 


t11)2(n,X,w)  = -^3  e X S wN(V)wN(WcosXtncosa,tn+sR^n>R(W 


n=l  s=l 


J oo  oo 

= (P~)2  f f wnKu)  wn(v)  cos  Xu  cos  wv  R(u)  R(y) 


N-l  N-n 


* /L  2 fn,n+s(u’’ 


(u,v)  du  dv 


n=l  s=l 


T j i 3(N,X,w)  = —p.  J J wN(u)  wN(v)  cos  Xu  cos  uv  R(u)  R(v) 


N-l  N-n 


•I  2 fn+s,n(u’y>du  dy 


n=l  s=l 


Now  (with  a change  of  variables  in  Tj  j j of  u=v  and  v=u) 

lTl)l)2(N>X.«)  + Tl,l,3(N,X,«)l 


OO  oo 

■ / / wN(u)  wN(v)  R(u)  R(v)  (cos  Xu  cos  cov  + cos  Xv  cos  ouu) 

A II 


N-l  N-n 


•2  2 fs(v-u)  fn(u)  dv  du 


n=  1 s=  1 


tv — 1 J 

OO  » ->  OO  rp 

J | R(u)  | 2,  fn(u)  / 1 R(y)  1 A fn(y-u)  dy  du 


oo 

0 J | R(u)  I du 
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We,  therefore,  have  the  following  lemma: 


Lemma  4.  Let  R(t)  eLj  n L2.  Then 


(00  \ 2 00 

0 I |R(u)|dul  + f R2(u)  du  + 0(\)  <f>(cj) . 
o ' o 

The  proof  of  Lemma  4 follows  (18)  through  (20). 

The  final  term  that  we  investigate  is  T^N.X.cj).  Before  stating  the  derivation  of  the 
asymptotic  behavoir  of  T4,  the  following  assumption  on  the  fourth  order  moment  of  the 
output,  hardlimited  signal  is  made. 


Assumption  1:  The  fourth  order  moment  function  of  the  output  signal, 


= e{y(0  Y(t+tj)  Y(t+t2)  Y(t+t3)}  , 


satisfies 


MY(tl>t2>t3)  I ^ n n b(tj) 


where  bj(t)  eLj  and  b(t)  eLj  and  bj,b  are  continuous,  even,  nonnegative  functions  which  are 
nonincreasing  over  [0,°°).  In  addition,  the  function  b(t)  satisfies  the  inequality 

b(tj)  •b(t2)<b(t1  + t2)  (21) 

for  tj  and  t2  positive. 

We  first  note  that  the  tj’s  in  the  above  assumption  are  real  numbers  and  not  the  random 
variables  from  the  Poisson  point  process  T.  Secondly,  the  inequality  (and  other  conditions) 
for  the  function  b(t)  in  (21)  is  satisfied  for  the  exponential  function  exp  (-a|tl). 


Lemma  5.  LetRy(t)eLj.  Then 

0(1+0)  I , J 

|T4(N,X,w)| max  |l,b(0)| 


The  proof  of  Lemma  5 comprises  Appendix  A. 


b(u)du j2  + (/  IRY<u>|duj2  +0(1) 


We  note  that  by  the  above  analysis,  i.e.,  by  combining  the  results  of  Lemmata  2 
through  5,  the  estimate  (5)  is  not  a consistent  estimate  of  the  true  spectral  density  ^(X). 
In  fact,  the  covariance  of  the  estimate  is  (asymptotically)  bounded  by 


/ oo  \ *•  oo 

CovJ?N(X),'0N(<o)|  = O fa  J I R(u) | du  j + J R‘ 

\ o / o 


(u)  du 


The  objective  of  future  work  is  to  obtain  tighter  bounds  and  show  that  these  new  bounds 
go  to  zero  asymptotically  as  N tends  to  infinity.  Specifically,  tighter  asymptotic  bounds 
need  to  be  found  for  the  two  terms  T|(N,X,<o)  and  T4  j(N,X,oj). 


■ 


CONCLUSIONS 

The  problem  of  reconstructing  the  power  spectral  density  of  a Gaussian  signal 
process  from  hardlimited  observations  taken  at  Poisson  sampling  instants  was  considered. 
The  estimate  of  the  power  spectral  density  of  the  original  (not  hardlimited)  signal  process 
is  shown  to  be  asymptotically  unbiased  as  the  number  of  observation  points  approaches 
infinity.  The  important  consequences  of  the  asymptotically  unbiased  estimator  are  as 
follows. 


1 . No  aliasing  of  the  spectral  density  of  the  (original)  signal  process  occurs  in  spite 
of  hardlimiting  the  amplitude  of  the  signal. 

2.  The  asymptotic  rate  of  convergence  of  the  bias  is  identical  to  the  rates  of  con- 
vergence for  the  cases  of  (regular)  Poisson  spectral  estimation  and  of  periodic  spectral 
estimation  (that  is,  the  convergence  of  the  periodic  estimate  to  the  aliased  spectral  density). 
In  particular,  under  mild  integrability  conditions  on  the  power  spectral  density,  the  rate  of 
convergence,  for  the  bias,  will  be  0(1 /N)  where  N is  the  number  of  observations. 

The  covariance  of  the  estimate  was  shown,  however,  only  to  be  bounded  as  the 
number  of  observations  tends  to  infinity.  Therefore,  the  crucial  task  that  remains  to  be 
completed  is  to  show  that  variance  of  the  estimate  tends  asymptotically  to  zero.  Upon  the 
completion  of  bounding  the  variance  of  the  estimate  by  a factor  which  tends  to  zero  as  the 
number  of  observations  tends  to  infinity,  the  estimate  of  the  spectral  density  will  then  be 
known  to  be  mean-square  consistent  — a fact  that  will  increase  the  confidence  of  the  esti- 
mate. This  task,  the  bounding  of  the  variance,  will  be  the  focus  of  future  work. 

In  addition  to  completing  the  proof  of  showing  that  the  estimate  is  mean-square 
consistent,  future  work  should  include  simulation  results,  for  finite  sample  size  N,  that 
would  indicate  the  feasibility  of  implementing  Poisson  spectral  estimation  using  hardlimited 
samples.  Furthermore,  since  each  hardlimited  sample  requires  only  one  bit  of  information 
for  the  purpose  of  storage  or  transmission,  a trade-off  study  between  the  decrease  in  the 
required  representation  (of  one  bit  per  sample)  which  decreases  storage  area  and  the 
(possible)  increase  in  the  number  of  samples  for  good  spectral  estimates  should  be 
conducted. 


3 

j 


l 
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APPENDIX  A 

Lemma  5 (restated  here  for  the  convenience  of  the  reader)  and  its  proof  comprise  the  infor- 
mation in  this  appendix. 


Lemma  5.  LetRy(t)eLj.  Then 
.0(1+0) 


|T4(N,\,cc)I<:  , 
402 


max  |l,  b(0)|^/  b(u)  duj  +(^j  |Ry(u)|duj 


+ 0(1). 


Proof.  To  begin  the  analysis  of  T4  (N,X,u>),  the  expected  value  with  respect  to  the  Y process 
will  be  taken  first  - only  the  relevant  terms  (with  respect  to  Y)  appear. 


R <N)(n)  - Ry(tn)  ] [ R^N>(k)  - Ry(tk)  | 

N-n  N-k 

= ^T  2 2 MY(tj  ~ ‘fi+n  “ lj+k  “ *£> 


8=1  j=l 

N-k 


- Ry(tn)/N  • ^ Ry^C+k  “ 4fi) 
8=1 


N-n 

2 
t=  i 


- Ry(tk)/N  • 2 RY(t8+n  ~ ^ 


+ Ry(tn)  Ry(tk) 
and  T4(N,\,w)  becomes 

T4(N,X,w)  = T4  , (N,X,w)  + T4j2(N,X,w)  I /402 


where 


T4>j(N,X,cu)  - Ej 


ANBNwN(tn)  wN(tk)c°s  Xtn  cos  cotk 


RY(tC+k  “ 4g) 


T4  2(N»  X,  w)  — E-p 


aNBN  wN^n^  wN^k^  cos  ^ *n  cos  wtk 


RY(tn)RY(tk)-RY(tk)/N-  I RY(tk+n-t8> 


The  term  T4  2(N,X,to)  is  now  considered  and  becomes 


f f WN^  WN^  cos 


T42(N,X,a;) 


COS  CJV 


Ry(u)  Ry(v)  - Ry(v)/N  • J Ry(u) 


fn  k(u,v)  du  dv 


wj^Cu)  wj^v)  cos  Xu  cos  cov 


7 1/2  7 1/2 

1 1 - R^(u)]  [ 1 - R2(v)]  Ry(u)  Ry(v)  n fn>k(u,v)  du  dv 


We  recall  from  the  proof  of  Lemma  2 that 


N N 


2 S fn,k<u>v 


n=l  k=l 


converges  as  N tends  to  infinity,  and  hence,  by  the  Kronecker  lemma  [61 


N N 


TT  S n S fR,k('i,v)^OasN 


n = l k=l 


and  by  the  Lebesgue  dominated  convergence  theorem 


N N 


lT4,2(N,X,<o)|  < J j Ry(u)  Ry(v)  X n S fn,k(u>v> du  dv 


n=l  k=l 


0 as  N -*■  °°  . 


We  now  consider  j(N,X,o>).  The  sum  over  the  indices  C and  j is  divided  into 
three  regions 

Rj=Jj £<j} 

R2={j,C:  j<fi<j  + k) 

R3={j,C:  j + k < fij  . 

In  R j , tj+jc  - tg  > tj+k  - tj  > 0 and,  therefore, 

b(tj+k  - tg)  < b(tj+k  - tj)  . (A- la) 

Also  note,  that  in  Rj 

b(tj  - tg)  < b(0)  . (A- lb) 

In  the  region  R2,  by  Assumption  1 

b(tg  - tj)  b(tJ+k  - tg)  < b(tj+k  - tj)  . (A-2) 


T 


*•  ; 


i- 


Finally,  in  the  region  R3,  tg-tj  > tJ+k  -tj  > <3,  and  the  following  inequalities  hold 

b(tj  - tg)  < b(tj+k  - tj)  (A.3a) 

b^j+k  “ < b(°)  • (A-3b) 


By  Assumption  1 and  the  inequalities  (A-l)  through  (A-3)  T4>  j (N,\,co)  becomes  bounded  by 

N N 


IT4,|(N,x,«)l  < 2 2 / / 


n=l  k=l  o o 


wN(u)  wn(v) 


N' 


ZZ  b(0)  b(u)  b(v)  + 2 2 b(u)  b(v)  + 2 1 b(0)  b(u)  b(v) 


R 


R~ 


Ri 


|RY(u)RY(v)l  (l-£) 


fn  k(u>v) du  dv  • 


(A-4) 


At  this  point  we  assume  a specific  covariance  averaging  kernel  is  given  - namely  the 
Dirichlet  spectral  window,  which  has  as  the  corresponding  covariance  averaging  kernel 


1 |tl<MN 

0 otherwise 


wN(t)  = h(t/MN)  = 

Thus,  (A-4)  is  bounded  by 

|T4>1(N,X,«)I<  2 2 / / 


mn  mn 


n=l  k = l o o 


*N-n  N-k  ‘ 

1 

w 

2 2. 

• 

.6=1  j = l . 

max  t l,b(0)]  b(u)b(v) 


+ |RY(u)RY(v)|  (l  -£■) 

mn  mn 

r-»  t-'  OO  OO 

- I l f f 


n=l  k=l  o o 


|RY(u)RY(v)|  (l-£) 


fn,k(u>v)  du  dv 
N2-(n+k)N  + nk 


N" 


max  [ 1 ,b(0)J  b(u)b(v) 


fn,k(u>v)  du  dv 


I 

1 1 

, I 
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Now,  since  ZZ  fn,k^u,v^  ^00  + 0)>  by  the  Lebesgue  dominated  convergence  theorem 


MN  mn 


OO  OO  r-s 

J J max  [ 1 , b(0)]  b(u)  b(v)  2,  A fn,k^u>v) du  dv 


o o 


n = l k=l 


<0(1  + 0)  max  [ l,b(0)] 


(OO 

1 


b(u)  du  | . 


J 


MN  mn 


OO  OO 

j J max  [ 1 , b(0)]  b(u)  b(v)  2,  2,  ° fn,k^u>v^ du  dv 


n=l  k=l 


<0(1  +0)  max  [ 1,  b(0)J 


>(/ 


b(u)  du 


-«(*) 


as  N -*■  «> 


For  the  third  term 


mn  mn 


OO  OO  « s rm 

[‘  J max  [ 1 , b(0)  ] b(u)  b(v)  2,  £ ^ ^(u.v)  du  dv 

n i N 


<0(1+0)  max  ( 1,  b(0)  ] 


n=l  k=l 


W (/ 


b(u)  du 


Finally, 


MN  mN 


/ J RY(u)  ry(v>  2 2 f1  "fi)  fn,k(u>v>dudv 

o o n=l  k=l  V ' 


<0(1+0)  If  Ry(u)  du  ] 


and  the  lemma  follows.  Q.E.D. 


(A-6) 


